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I. INTRODUCTION 



When a quantum field theory is constructed on a curved background spacetime, the gravitationally 
induced vacuum polarization typically induces a non-zero vacuum expectation value for the quantum 
stress-energy tensor. It is well-known that this gravitational vacuum polarization often violates 
the point-wise energy conditions and sometimes violates the averaged energy conditions. In this 
paper I begin a systematic attack on the question of where and by how much the various energy 
conditions are violated. To keep the discussion manageable, I work in the test-field limit, and focus 
on conformally coupled massless scalar fields in Schwarzschild spacetime, using the Hartle-Hawking 
vacuum. The discussion invokes a mixture of analytical and numerical techniques, and critically 
compares the qualitative behaviour to be expected from the Page approximation with that adduced 
^ from the numerical calculations of Anderson, Hiscock, and Samuel. I show that the various point- 

^ I wise energy conditions are violated in a series of onion-like layers located between the unstable 

, photon orbit and the event horizon, the sequence of violations being DEC, WEC, and (NEC+SEC). 

' Furthermore the ANEC is violated for some of the null geodesies trapped in this region. Having 

CO ' established the basic machinery in this paper, the Boulware vacuum will be treated in a companion 

paper, while studies of the Unruh vacuum should be straightforward, as should extensions to non- 
conformal couplings, massive scalars, and Reissner-Nordstrom geometries. 

> 

o 
o 
^ 
o 

' When a quantum field theory is constructed on a curved background spacetime, the gravitationally induced vacuum 
polarization typically induces a non-zero vacuum expectation value for the stress-energy tensor [yj-|6|. It should be 
^ • pointed out that fully self consistent calculations are horrendously difficult; consequently most calculations in the 
literature are carried out in the test-field limit. (Wherein the vacuum polarization is not permitted to back-react 
r* on the geometry via the Einstein field equations.) Even in the test-field limit, most known results are obtained by 
. !^ I numerical (rather than analytic) computations. 

. Perturbative self-consistent solutions around flat spacetime have recently been investigated by Flanagan and 
^ ■ Wald , building on earlier work of Horowitz |^ , and Horowitz and Wald |^ . 

Another type of energy condition, based on the "quantum inequalities" , is investigated for a Schwarzschild back- 
ground in recent papers by Ford and Roman [ p^ , pd| . 

Ford and Roman have also discussed the ANEC and AWEC in (l+l)-dimensional and (3-|-l)-dimensional evapo- 
rating black holes (Unruh vacuum) pl , and for the (1-fl) -dimensional Boulware vacuum. In contrast, the present 
paper works strictly in (3-1-1) dimensions and treats the equilibrium Hartle-Hawking vacuum state. 

Now we know, on rather general grounds (the existence of the Hawking radiation effect and concomitant violation 
of the classical area increase theorem of black hole dynamics) that at least some of the energy conditions must be 
violated at or near the event horizon of a black hole. In this paper I shall explore these issues in a little more detail. 
I shall continue to work in the test-field limit: I discuss the effect of the gravitational vacuum polarization on various 
energy conditions. In particular, I discuss the point-wise null, weak, strong, and dominant energy conditions (NEC, 
WEC, SEC, and DEC), the averaged null energy condition (ANEC), and shall furthermore introduce and discuss 
several new, potentially interesting, energy conditions (PNEC, Scri-NEC, and Scri-ANEC). 

For the purposes of this paper, I restrict attention to the most well-studied curved-space quantum field theory: 
the conformally-coupled massless scalar field on Schwarzschild spacetime, in the Hartle-Hawking vacuum. For this 
geometry and vacuum state one has both (1) a useful analytic approximation to the gravitational polarization — 
Page's approximation |l2| , and (2) numerical estimates of the vacuum polarization — see the numerical calculations of 
Howard [Of, Howard and Candelas ||l3], and Anderson, Hiscock, and Samuel |]l5|-0|.) 
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Having established the basic machinery in this paper, a study of the Boulware vacuum will be presented in a 
companion paper, while further studies of the Unruh vacua should be straightforward, although tedious. (For the 
current state of affairs in the Unruh vacuum, see Similarly, extensions to non-conformal couplings, massive 

scalars, and Reissner-Nordstrom geometries should be straightforward. 



II. NULL ENERGY CONDITIONS 



To set the stage, recall some basic definitions — 
Definition [NEC] 

The null energy condition is said to hold at a point p if for all null vectors 



T^, k'^k" > 0. (1) 



Definition [ANEC] 

The averaged null energy condition is said to hold on a null curve 7 if 

T^^ k^k" dX > 0. (2) 

Here A is a generalized affine parameter (see [^8[ pp. 259, 278, 291]) for the null curve, the tangent vector being 
denoted by fc'^. 

Note that in almost all applications it suffices to consider null geodesies, rather than gen eric null curves. However, 
isolated discussions of the ANEC on non-geodesic null curves do occur in the literature — the extended definition 
presented here allows one to compare and contrast the present results with other calculations. 

As a practical matter I shall often replace the word "inextendible" by the phrase "inextendible past the event 
horizon" . This is a purely pragmatic decision based on a lack of numerical data inside the event horizon, coupled 
with the general feeling in the community that Page's analytic approximation will become progressively worse as one 
approaches the central singularity. (Thus my studies of the ANEC can more precisely be referred to as studies of a 
"Truncated- ANEC".) 

In an earlier publication |20|, (see also j|]), I derived a no-go theorem for the ANEC: I showed that the ANEC 
is violated (in the test-field limit) whenever the background spacetime has a non-zero scale anomaly. Specifically, I 
proved the 

ANEC no-go theorem 

In (3+ 1) -dimensional spacetime, {M,g), for any conformally coupled massless quantum field, in any conformal 
quantum state: If 

C'^p. ^ (3) 

then, holding the renormalization scale /i fixed, it is possible to find a rescaled metric g — fl^g such that the ANEC is 
violated on the spacetime (M,g). 

This no-go theorem, since it depends on a special case of the quantum field theoretic (conformal) anomaly, is known 
to be independent of the choice of vacuum state. (One can easily generalize this argument to include non-conformally 
coupled fields, see 0.) 

Unfortunately, a brief calculation suffices to show that the tensor vanishes on Schwarzschild spacetime. Indeed 
it vanishes on any Ricci-flat spacetime (or indeed on any spacetime conformal to an Einstein spacetime). Thus 
Schwarzschild spacetime avoids the no-go theorem enunciated above, and it becomes a matter of explicit checking 
to see whether or not the ANEC is satisfied or violated on this geometry. Note that in doing this explicit checking 
we will have to be cognizant of our particular choice of vacuum state, and indeed our results will also depend on the 
particular class of null geodesies under consideration. 

In performing the explicit checks alluded to above, I have found it useful to define several new energy conditions 
that are in some sense intermediate between the NEC and ANEC: 
Definition [PNEC] 

The partial null energy condition will be said to hold at a point p for a set of null tangent vectors P{p) C Tp if for all 
null vectors k^ G P{p) 

T^, k^k-' > 0. (4) 
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Now, if one does not place any further constraints on the set the resulting definition, while general, is not 

particularly useful. (It is best used as a diagnostic tool to probe interesting regions of the tangent bundle.) A suitable 
restriction is: 
Definition [Scri NEC] 

The asymptotic null energy condition will be said to hold on an asymptotically flat spacetime if for every point p of 
the spacetime, and all null tangent vectors such that the associated null geodesic through p either (1) escapes to 
future null infinity (Scri+), or (2) arrives from past null infinity (Scri~), one has: 

T^, k'^k'' > 0. (5) 



(This definition basically says that one should not worry too much if violations of the NEC occur only when one 
looks along null geodesies that never make it out to null infinity.) 

We shall soon see that in the Hartle-Hawking vacuum the gravitational vacuum polarization of a conformally- 
coupled massless scalar field residing on the Schwarzschild background geometry everywhere satisfies this Scri~NEC 
energy condition, though it does not everywhere satisfy the NEC. We shall also use the PNEC to investigate the 
regions in phase space (the tangent bundle) where suitable positivity conditions hold. 

Finally, for completeness, I enunciate 
Definition [Scri ANEC] 

The asymptotic averaged null energy condition will be said to hold on an asymptotically flat spacetime if every in- 
extendible null curve which either (1) escapes to future null infinity (Scri+), or (2) arrives from past null infinity 
(Scri"), satisfies the ANEC: 

T^^ k^'k" dX > 0. (6) 

(This definition basically says that one should not worry too much if violations of the ANEC occur only along 
null geodesies that never make it out to null infinity. After all, null infinity is where you want to work to establish 
results such as the Penrose-Sorkin-Woolgar positive mass theorem or the Fricdman-Schleich-Witt topological 
censorship theorem p2|.) 

From the previous comment regarding Scri-NEC it automatically follows that, in the Hartle-Hawking vacuum, 
the gravitational vacuum polarization of a conformally coupled scalar field residing on the Schwarzschild background 
geometry satisfies this Scri-ANEC. We shall also see that not all null geodesies satisfy the ANEC itself. 



III. VACUUM POLARIZATION IN SCHWARZSCHILD SPACETIME: HARTLE-HAWKING VACUUM 



In the particular case of conformally coupled scalar fields residing on Schwarzschild spacetime one has a lot of 
information about the vacuum polarization. By spherical symmetry one knows that 
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Where p, r and p are functions of r. Ad and h. Note that I set G = 1, and choose to work in a local-Lorentz basis 
attached to the fiducial static observers (FIDOS). 

Page's analytic approximation gives |12[l^ a polynomial approximation to the stress-energy tensor: 



p{r) = +3po 
T(r) = -poo 
p{r) = +p,3o 
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Here I have defined a constant (to be interpreted as the pressure at spatial infinity) by 



Po 



90(167r)2(2M)^ 



(12) 



Note that I have exphcitly expanded the functions given in references as polynomials in 2M/r to exhibit the 

fact that all components of the stress-energy are well behaved at the horizon. It is perhaps somewhat surprising 
that the rather messy rational quotients of polynomials exhibited in [p"2|jr^ ] reduce to such relatively nice compact 
expressions. These formulae have also been checked for consistency with the Brown-Ottcwill extensions of the original 
Page approximation see esp. p 2517]. See also Elster If one prefers, an alternative form for the Page 

approximation is 



{H\T^o\H)=Poo {P{r) 



The polynomials P{r) and Q{r) being 
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and 



^, , ^ ^f2M\ f2My f2MY f2MY f 2My f2M\ 



Q{r) = 24 



\ r J 



The trace of the stress-energy tensor is given by 



{H\T^f,\H) = 9&p, 



(14) 



(15) 



(16) 



This result, because it is simply a restatement of the conformal anomaly, is known to be exact. 



At the event horizon, Page's analytic approximation gives |12 



{H\T^o\H)=p„ 
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The explicit numerical calculations of Howard [^, Howard and Candelas [|^, and Anderson, Hiscock, and 
Samuel |p^-^ show that Page's analytic approximation is reasonably good — the worst deviations occur in the im- 
mediate vicinity of the event horizon where the Anderson-Hiscock-Samuel numerical analysis gives 



{H\T^o\H) 
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(18) 



A convenient table of numerical results, covering the range r — 2M to r — 6.8M may be found on page 2541 of 
reference [|l^. (Near the event horizon, the Howard-Candelas data is believed to be accurate to two significant 
digits. The accuracy is expected to improve with increasing r, as the numerical results converge upon the Page 
approximation.) 

Additionally, I have obtained via private communication the more intensive numerical data of Anderson-Hiscock- 
Samuel [p5|-p7[, which covers the range r = 2M to r = 5M at much higher resolution. (Near the event horizon, the 
Anderson-Hiscock-Samuel data is believed to be accurate to three significant digits.) I have collated the data and 
used it to construct suitable interpolating functions to be used for the numerical aspects of the following discussion. 

To construct the interpolating functions I use the Anderson-Hiscock-Samuel data in the range r = 2M to r — 5M, 
and the Howard-Candelas data in the range r = 5M to r = 6.8M, augmented by the fact that we know the exact 
result at r = oo. Using Mathematica, I have fit the data using third-order interpolating polynomials in the variable 
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z = 2M/r. (Note that z conveniently runs from z = 1 at the horizon to z = at spatial infinity; z is also the natural 
variable to use when numerically analyzing the Page approximation.) 

Of course, there are may other pape rs extant which calculate the stress-energy tensor in somewhat different con- 
figurations. See, for example, [n2 



15 



1( 



|,|23|-|34| . More recently, in a very interesting development, Matyjasek |3^] has 
developed a curve-fitting analysis that fits the numerical data to high accuracy. In this paper I prefer to work directly 
with the numerical data. In this paper I am only presenting the minimum requirements for the particular job at hand 
and it is clear that with additional effort more information can be extracted by mining this additional vein of data. 



IV. NULL ENERGY CONDITION 



A. Outside the horizon: 



Outside the event horizon, the NEC reduces to the pair of constraints 

p{r) - r(r) > 0? p{r) + p{r) > 0? 

Page's approximation yields 



(19) 



p{r) ~ T{r) = 4pc 
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(Note the factorization!) Consequently p(r) — T(r) is always explicitly positive outside the horizon. Furthermore 



p{r) ^ p{r) = Apo 
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(21) 



Numerically finding the roots of this sixth order polynomial shows that p(r) + p(r) is negative from r = 2M to 
r « 2.18994Af . 

Thus the Page approximation suggests that the null energy condition is violated in the range r G [2M, 2.18994M]. 

Turning to the numerical data, one easily verifies that that outside the horizon p{r) — T(r) > 0, while p = r at the 
horizon. This is in complete agreement with Page's analytic approximation. 

On the other hand, inspection of the numerical data indicates that p -l-p < for r < 2.298 A/. (This number was 
obtained by fitting third-order interpolating polynomials to the numeric data; and then numerically finding the root.) 
This is qualitatively, though not quantitatively in agreement with Page's analytic approximation. 

Thus the numerical data indicate that the null energy condition is violated in the range r € [2M, 2.298M]. 



B. Inside the horizon: 



Inside the event horizon, the radial coordinate becomes timelike, and the roles played by p(r) and T(r) are inter- 
changed. The NEC reduces to the pair of constraints 

T(r) - p(r) > 0? T(r) + p{r) > 0? (22) 

Unfortunately, inside the horizon I am reduced to reliance upon the Page approximation by a total lack of numerical 
data. Now there are some subtleties associated with extrapolating the Page approximation inside the event horizon. 
The existence of a Killing vector that points in the t direction is absolutely critical, (even if the t direction is no longer 
timelike). Thus if you wish to extrapolate the Page approximation inside the event horizon it is necessary to limit 
attention to an eternal black hole (that is, the maximally extended Kruskal-Szekeres manifold). 

I should remind the reader that in the Hartle-Hawking vacuum the stress-energy is at least known to be regular at 
the horizon, so the Page approximation should not be too far wrong just inside the event horizon of an eternal black 
hole. On the other hand, one does expect the Page approximation to get progressively worse as one moves further in 
toward the singularity, so one should probably not take these results too seriously once one is far inside the horizon. 

I should especially warn the reader against extrapolating the Page approximation to the interior of an astrophysical 
black hole: we really do not expect the interior of an astrophysical black hole to have the same Killing vectors as the 
exterior, but instead expect the interior to be fully dynamical. (In addition we do not expect an astrophysical black 
hole to even to be in the Hartle-Hawking quantum state.) 
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The present calculations inside the horizon are strictly limited to eternal black holes and are exhibited merely 
because they arc do-able, and because they may be suggestive of the actual situation. 

Inspection of Page's analytic approximation indicates that T(r) — p{r) > for all r < 2M, so this part of the NEC 
is satisfied inside the event horizon. On the other hand 

r{r)+p{r) = -2Ap^( — \ , (23) 

which is everywhere negative, (both inside and outside the horizon). Thus the Page approximation suggests that the 
null energy condition is violated throughout the interior of the black hole. 

This is enough to inform us that, (assuming the reliability of the Page approximation in this matter), the weak, 
strong, and dominant energy conditions are also violated throughout the entire interior of the black hole. 

V. WEAK ENERGY CONDITION 
A. Outside the horizon: 

Outside the horizon, the weak energy condition is equivalent to the three constraints 

pir) > 0? 

p{r) - r(r) > 0? 

p(r)+p(r)>0? (24) 

The last two constraints have already been discussed — they are simply the null energy condition. In the Page 
approximation, finding the roots of the relevant sixth-order polynomial shows that the density constraint is violated 
for r < 2.3468M. Therefore the weak energy condition is violated in the range r G [2M, 2.3468M]. 

In contrast, the numerical data indicate that violations of the weak energy condition extend to the region r G 
[2M,2.438M]. 

B. Inside the horizon: 

Inside the horizon, the weak energy condition is equivalent to the three constraints 

r(r) > 0? 
r(r) - p{r) > 0? 

T{r)+p{r)> 07 (25) 

The last two constraints have already been discussed (when dealing with the NEC). The violation of the NEC is 
already sufficient to tell us that the WEC is violated everywhere inside the horizon, but we can add in passing that 

the Page approximation additionally implies that T(r) is everywhere negative inside (and outside) the horizon. 

The Page approximation suggests that the weak energy condition is violated throughout the interior of the black 
hole. 

VI. STRONG ENERGY CONDITION 
A. Outside the horizon: 

Outside the horizon, the strong energy condition is equivalent to the three constraints 

p{r) - T(r) > 0? 
p{r) +p{r) > 0? 

p{r) - T(r) + 2p{r) > 0? (26) 

We have already looked at the first two constraints when discussing the null energy condition. The third condition is 

always satisfied outside the horizon (since both p{r) — T(r) > 0, and p(r) > in this region). This is true both in the 
Page approximation, and by appeal to the numerical data. Thus the strong energy condition is violated in the same 
region as the null energy condition. 
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B. Inside the horizon: 



Inside the horizon, the strong energy condition is equivalent to the three constraints 

T(r) - p{r) > 0? 
T(r) +p{r) > 0? 

T(r) - p{r) + 2p{r) > 0? (27) 

We have aheady looked at the first two constraints when discussing the null energy condition, and thereby know that 
(the Page approximation suggests that) the strong energy condition is violated throughout the entire interior of the 
black hole. 

For completeness I point out that the first condition is satisfied throughout the interior, the second of these conditions 
is violated throughout the interior, while the third condition also fails throughout the interior. 



VII. DOMINANT ENERGY CONDITION 

A. Outside the horizon: 

Outside the horizon, the dominant energy condition is equivalent to the three constraints 

p(r) > 0? 

r(r) e [-p(r), +p(r)]? 

p{r)€[-p{r),+p{r)]? (28) 

We can rephrase this as 

p{r) > 0? p{r) ± T(r) > 0? p{r) ± p{r) > 0? (29) 
Using the Page approximation, and restricting attention to the range [2M, oo] , one has: 

p(r)<0 r e [2M,2.3468M]; 

p{r) - T{r) < r G 0; 

p{r) + T(r) < re [2M, 2.77256M]; 

p{r) - p{r) < r G [2M, 2.58512M]; 

p{r) +p{r) < r G [2M, 2.18994Af]. (30) 

Pulling this all together, the Page approximation suggests that dominant energy condition fails in the region r G 
[2M,2.77256M]. 

The numerical data implies quantitative though not qualitative modifications. In the range [2M, oo] one finds: 

p(r)<0 r G [2M,2.438M]; 

p{r) - T{r) < r G 0; 

p{r) + r(r) < r G [2M, 2.992M]; 

p{r) - p{r) < re [2M, 2.628M] ; 

p{r) +p{r) < re [2M,2.298M]. (31) 

Pulling this all together, the numerical data indicate that the dominant energy condition fails in the region r e 
[2M, 2.992M]. Note that this is suspiciously close to r = 3M — the unstable circular photon orbit — and might be 
trying to tell us something. 

If one actually calculates {p+t) / {\p\-\-\t\) at r = 3Af one gets 6.02x 10^'^ given the expected three-significant-digit 
numerical reliability of the data this is (questionably) compatible with zero. If one takes this suggestion seriously, it 
would imply a hidden (accidental?) symmetry in the stress tensor at r = 3M: 



{H\T~^'\H)l^,M^ 
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6 



(32) 
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It is quite possible however, that this is purely a numerical accident. 

Finally, I point out that this conjectured property certainly does not survive the introduction of non-conformal 
coupling, nor is there any particular reason to expect it to survive the introduction of non-zero rest mass. This 
conjecture also most definitely does not hold in the Boulware or Unruh vacuum states ||3^ , |3^ , |lT| . 



B. Inside the horizon: 



Inside the horizon, the dominant energy condition is equivalent to the three constraints 

r(r) > 0? 

p(r) e [-T(r),+T(r)]? 

p(r) e [-T(r),+T(r)]? (33) 

We can rephrase this as 

-r(r)>0? -r±p(r)>0? - r ± p(r) > 0? (34) 

We already know that r(r) < inside the horizon (in fact everywhere). So, provided the Page approximation is not 
misleading in this regard, the dominant energy condition is violated throughout the interior of the black hole. No 
additional information comes from the other conditions. 



VIII. PARTIAL NULL ENERGY CONDITION 



A. Outside the horizon: 



To analyse the partial null energy condition introduced earlier in this paper, consider a generic null vector inclined 
at an angle ip away from the radial direction. Then without loss of generality, in an orthonormal frame attached to 
the {t, r, 0, (p) coordinate system, 

/c'^ cx (±l,cos?A,0,sin?A). (35) 

Ignoring the (presently irrelevant) overall normalization of the null vector, one has 

Ti.^k^k" (x{p-T cos^ ip+p sin^ ^p) 

= i[p-T] + [T +p]sin^^P). (36) 

We have already seen that p — r is positive outside the event horizon. On the other hand the Page approximation 
gives 

r(r)+p(r) = -24poo f— j , (37) 

which is everywhere negative. A glance at Howard's numerical data confirms that the numerical data also satisfies 
T +p <0. 

Thus the partial null energy condition fails for those radii r, and those angles V', such that 



sini; > sin[V^,Ht(r)] ^ J ^fl-^. (38) 
V \r{r) +p{r)\ 

That is, the partial null energy condition fails for ijj e [V'crit(?'), ■'i'/2]. Defining z = 2M/r this critical angle is (in the 
Page approximation) given by 

ipcritiz) = sin^i - + 3z + 6z2 + lOz^ + 15z4 + 21z^)/6z^] . (39) 

For r > 2.18994M {z < 0.913267) there are no real solutions to this equation. At r = 2.18994M one has Vcrit = 7r/2, 
while V'crit('") moves monotonically to zero as r ^ 2M. (See figure |l| where ipcrit{z) is plotted as a function of z.) 
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Thus, sufficiently near the horizon, almost all directions violate the partial null energy condition. As one moves 
further away from the horizon the violations of the partial energy condition are confined to null vectors that are 
progressively more and more transverse, finally at r = 2.18994Af the violations of the partial null energy condition 
disappear completely. 

Using the numerical data will modify the precise location where this behaviour manifests itself, but will not qualita- 
tively modify this picture. The violations of the partial null energy condition vanish outside r ~ 2.298Af {z — 0.8703), 
and the numerically determined values of ipcntiz) are superimposed on figure |^. 




0.86 0.88 0.9 0.92 0.94 0.96 0.98 



FIG. 1. PNEC: Outside the horizon. This graph shows the (normalized) critical angle ■!/;crit(^)/[7r/2] above which the partial 
null energy condition fails. (The solid line represents the numerical data; the dotted line represents the Page approximation.) 



B. Inside the horizon: 



Inside the event horizon, one should consider a generic null vector inclined at an angle V' away from the t direction 
(which is now spacelike). Then without loss of generality, in an orthonormal frame attached to the (t, r, 9, cj)) coordinate 
system. 



cx (cos ■0, ±1, 0, sinT/i). 



One should now consider the quantity 



T^^k^k" (X {t — p cos^ ip + p sin^ -0) 
= {[t- P] + [P + P] sin^Vi). 



(40) 



(41) 



Note that inside the horizon r — p is strictly positive, while p + p is strictly negative. In fact, the Page approximation 
yields 



p(r) + p(r) = Apo 



f2M\ f2MY f2MY f2MY f2MY f2MY 

:y-)y-) +n-) ^i-) 

One deduces that violations of the partial null energy condition occur whenever 

sin?/; > sin[i/icrit(?')] ; 



(42) 



r(r) - p{r) 



\p{r)+p{r)\ 



(43) 



Consequently, violations of the partial null energy condition arc confined to the range € ['!/'crit(''), 7r/2]. Again 
defining z = 2M/r this critical angle is (in the Page approximation) given by 



V'crit(2) = sin 1 



' (1 - z)(l + 3z + 6z2 + 10z3 + Ibz* + 21^5) 
(1 + 2z + 3z2 + 4z3 + 5z4 + 6z5 - 27z^) 



(44) 
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One has, at r = 0, 'tjjcrit{r = 0) = sin~^ 1^7/9 « 62deg, while i^cmt falls to zero as r approaches 2M. Qualitatively: 
near (but inside) the horizon almost all null directions suffer violations of the partial energy condition, while near the 
central singularity somewhat fewer directions violate the partial null energy condition. 

Again, I remind the reader to not take the Page approximation too seriously as one approaches the singularity. 



IX. ASYMPTOTIC NULL ENERGY CONDITION 



I now turn attention to the asymptotic null energy condition (Scri-NEC) that I introduced earlier in this paper. 
Discussing this energy condition will require some standard results concerning the null geodesies of the Schwarzschild 
spacetime. (See, for instance, such standard textbooks as Misner-Thorne- Wheeler Wald |3^], or Weinberg p8[.) 

There are two conserved quantities for null geodesic motion in Schwarzschild spacetime; the energy and angular 
momentum. Using these conservation laws the affine parameter can be chosen in such a way that 



dt 
d\ 

dX 



1 



1 - 2M/r ' 

2 



a 



(45) 
(46) 



The parameter a is the angular momentum per unit energy. If the null geodesic reaches asymptotic spatial infinity 
then this parameter can also be interpreted as the "impact parameter" . However, there is a large class of null geodesies 
that never reaches spatial infinity, for these null geodesies the notion of "impact parameter" is at best an abuse of 
language. 

The angle between the null geodesic and the radial direction is given by 



We are ultimately interested in the quantity 

T^.fc^fc" qc{[p-t] + [t+p] sin2 tf,) 

^{[p-t] + [r + p]{l-2M/r)ayr^}, 

but will need to go through a few preliminaries. 
The radial motion of null geodesies is governed by 



(47) 



(48) 
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1 - (1 - 2M/r)- 



So the turning points, dr/dt — 0, are given by the cubic 

= 0^(1 - 2M/r). 



(49) 



(50) 



If the impact parameter is small, a < 3\^M, then it is a standard result that there are no turning points: the null 
geodesic either plunges into the future singularity, or emerges from the past singularity of the maximally extended 
Schwarzschild spacetime. (For a S SVSM the geodesic may make a large number of "orbits" before crossing the event 
horizon.) 

If the impact parameter is marginal, a — SVSM, then it is a standard result that r = 3M, corresponding to the 
(unstable) circular photon orbit. 

If the impact parameter is large, a > SVSAd, then it is a standard result that there are two turning points at 
physical values of r. One of these turning points lies in the range r S (3M, oo), while the other lies in the range 
r e (2M,3M). 

Note that if a ^ 2M then the three mathematical roots of the cubic are approximately r « ±a — M, and 
r « 2M[1 + (2M/a)2]. The two physical roots are approximately r « a - M and r « 2M[1 + (2M/a)^]. 

The first of these turning points, r € (3M, oo), corresponds to the obvious class of null geodesies with high impact 
parameter — those geodesies that come in from spatial infinity, bounce off the angular momentum barrier, and return 
to spatial infinity. (For a > 3^/3M the geodesic may make a large number of "orbits" before returning to spatial 
infinity. ) 
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The second of these turning points, r G (2M, 3M), corresponds to a completely separate class of null geodesies with 
high "impact parameter" — these geodesies emerge from the event horizon at t = — oo, with high angular momentum, 
make a large number of "orbits" before reaching their maximum height above the event horizon, and then make an 
equally large number of "orbits" before returning to re-cross the event horizon at i = +00. For these geodesies the 
use of the phrase "impact parameter" to describe the parameter a is most definitely an abuse of language. 

To analyze the asymptotic null energy condition, start by first considering the quantity 



n(r,a) = (p 



' ip + psin^ 



(51) 



which I shall refer to as the NEC density. (Remember that 'ip{r, a) is an explicitly known function of r and a.) 

Those null geodesies that come in from infinity and return to infinity never get closer to the origin than r — 3M. 
See, for instance, [^6| , pages 672-678]. Inspection of either the Page approximation or of the numerical data indicates 
that p, — T, and p are all positive for r > 3M. Thus the asymptotic null energy condition is satisfied along all null 
geodesies that come from, and return to, infinity. 

For other classes of null geodesies it proves convenient to rewrite the quantity 11 (r, a) as: 



n(r,a) = {[p-t] + [T+p]{l-2M/r)ayr^}. 



(52) 



Now, for incoming null geodesies with smaller than critical impact parameter, the null geodesic may circle the black 
hole a large number of times, but is guaranteed to ultimately plunge into the event horizon pages 672-678]. This 
makes the analysis a little more subtle. Inspection of either the Page approximation or the numerical data shows that 
(outside the horizon) p — r is always positive, while r + p is always negative. Now use the fact that for an infalling 
null geodesic a < 3\/3M. Since t + p is negative, this implies (for this class of null geodesies) a lower bound 

n(r, a) > L{r) = {[p - r] + [r + p](l - 2M / r)27 M"^ / r^} . (53) 

Inspection of Page's approximation indicates that this lower bound is strictly positive (zero at the horizon) . In fact 



L(r)=4poo(l-^ 



f2M\ /2mV /2mV /2mV /2mV 81 /2A/ 



(54) 



(Note the factorization!) Thus the Page approximation suggests that the asymptotic null energy condition holds on 
all infalling null geodesies. (At least until one crosses the horizon!) This lower bound L{r) is plotted in figure 

By time reversal, this suggests that the asymptotic null energy condition also holds on all outgoing null geodesies 
that reach infinity. 




FIG. 2. L(z): A bound on NEC violations. For all impact parameters a < SVSM, the NEC density Tl{r,a) is bounded 
below by the quantity L{r), which is itself bounded below by zero. This implies that Scri-NEC is satisfied on all infalling or 
outfalling null geodesies. (The solid line represents the numerical data; the dotted line represents the Page approximation.) 
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With respect to the numerical data, one can easily evaluate L{r) numerically. The results are superimposed on 
figure 1^. While there are marked differences between the analytic approximation and the numerical data, both curves 
are seen to be strictly positive outside the horizon. 

We conclude that Scri-NEC is satisfied. For all null geodesies that reach spatial infinity the NEC density n(r, a) is 
strictly positive everywhere outside the event horizon. 

For completeness, and future reference, I point out that the Page approximation yields 



2M 

n(r,a) =4poo ( 1 



/2M\ /2mV /2mV /2mV /2mV f 2M 



(55) 



Again, note the factorization! 



X. TRAPPED NULL GEODESICS: 



Now turn attention to the trapped null geodesies. (These are null geodesies with high "impact parameter" a > 
SVSM, trapped in the region r g [2M,3M].) While these trapped null geodesies are, by definition, not directly 
relevant to the Scri-NEC, the tools developed above permit us to gain additional insight into the PNEC on trapped 
null geodesies. 

Pick some value of r in the range (2M, 3M). Then 11 (r, a) is guaranteed to be negative if one chooses 



'{p-t) 



{l-2M/r)\T+p\ 



2M, 



(p-r) 



Z2(l. 



Z T- 



(56) 
(57) 



(Setting z — 2M/r is again a convenient choice of variables for both analytic and numeric work. In this section we will 
only be interested in the region z S [2/3, 1].) Note that what we are doing is guaranteeing that with these definitions 
the PNEC is violated for the region of the (a, z) plane above the curve acrit(-z)- 
Using Page's approximation, this critical impact parameter is given by 



acrit(z) = 2M^J{1 + 'iz + 6z2 + I0z3 + i^z^ + 21z5)/(6z8). 



(58) 



The critical impact parameter implied by the numerical data was also determined and both curves are plotted on 
figure pi 




0.7 0.75 0.8 0.85 0.9 0.95 
FIG. 3. acrixiz) / M : The critical impact parameter above which NEC is violated on trapped null geodesies. (The solid line 
represents the numerical data; the dotted line represents the Page approximation.) 
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It should be noted that, for a given value of r G (2M, 3M), one cannot choose a arbitrarily: In order for a trapped 
null geodesic with impact parameter a to ever reach radius r it is necessary that a be small enough. Indeed one must 
have 



a < a„ 



2M 



2M/r ^z-^(\-z) 



(59) 



This may be thought of as a kinematic bound on trapped null geodesies: the region of the (a, z) plane above the 
curve amax(-z) is kinematically inaccessible to trapped null geodesies. (Note that the region below a — SVSM is 
also kinematically inaccessible.) The relevant curve is plotted in figures^ and where it is overlain with acritiz) as 
obtained from figure |^. By looking at where these curves cross one another we can draw some general conclusions. 
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0.7 0.75 0.8 0.85 0.9 0.95 
FIG. 4. Page approximation: The critical impact parameter above which NEC is violated (acrit(z)/M)is superimposed 
on the kinematic bound (amax(z)/M). (The dotted line represents the Page approximation; the dashed line represents the 
kinematic bound.) 



In the Page approximation, now restricting attention to the region outside the event horizon: 

1. Some trapped null geodesies (those with a £ [3V3M,2y'Yj3M] = [5.19615M, 6.1101M]) never experience NEC 
density violations. That is: n(r, a) > is satisfied along the entire portion of the null geodesic that lies outside 
the event horizon. (And therefore the [truncated] ANEC must be satisfied on this entire class of null geodesies.) 

2. Some trapped null geodesies (those with a > 7.436M) always experience NEC density violations. That is: 
n(r, a) < along the entire [truncated] null geodesic. (And therefore the [truncated] ANEC must be violated 
on this entire class of null geodesies.) 

3. All other trapped null geodesies (those with a E [6.1101M, 7.436M]) will experience some NEC violations when 
they get sufficiently close to the event horizon. (And therefore investigating ANEC violations on this class of 
null geodesies requires more work.) 

4. Note that a = 6.1101M corresponds to a null geodesic that reaches a maximum radius rmax ~ 2.18994M, a 
number that we have seen before (when discussing the NEC and PNEC). 
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0.7 0.75 0.8 0.85 0.9 0.95 
FIG. 5. Page approximation: The critical impact parameter above which NEC is violated {acTit{z)/M)is superimposed on 
the kinematic bound (amax(2)/Af). (The solid line represents the numerical data; the dashed line represents the kinematic 
bound.) 

Use of the numerical data implies quantitative though not qualitative changes. We observe that 

1. The set of trapped null geodesies which never experience NEC density violations is much smaller — those with 
a G [3V3M, 5.276M] = [5.19615M, 5.276M]. (The [truncated] ANEC must be satisfied on this entire class of 
null geodesies.) 

2. Some trapped null geodesies (those with a > 6.383M) always experience NEC density violations. (Therefore 
the [truncated] ANEC must be violated on this entire class of null geodesies.) 

3. All other trapped null geodesies (those with a € [5.276M, 6.383Af]) will experience some NEC violations when 
they get sufficiently close to the event horizon. (And therefore investigating ANEC violations on this class of 
null geodesies requires more work.) 

4. Note that a — 6.383Af corresponds to null geodesic that reaches a maximum radius rmax = 2.298A/, a number 
that we have seen before (when discussing the PNEC). 

The message to be extracted is this: The numeric data and analytic approximations are in good qualitative agree- 
ment with each other. Null geodesies that reach null infinity are well behaved (Scri-NEC is satisfied), but some of 
the trapped null geodesies encounter NEC violations (and ANEC violations). 

XI. AVERAGED NULL ENERGY CONDITION 

We are finally in a position to pin down precisely violations of the ANEC itself. Although significant information 
can already be extracted by using the point-wise energy conditions already discussed, beyond a certain stage explicit 
use of the af&ne parameterization must be invoked. Write the ANEC integral as fi page 117] 



'7 

{p — T cos^ ip + p sin^ ip) dX, 
n(r, a) ^2 dX. 



(60) 
(61) 
(62) 
(63) 



Note that the integrand is proportional to the quantity n(r, a) which has already been extensively discussed in the 
context of the Scri-NEC. 
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A. Scri-ANEC 



The arguments previously adduced for the Scri-NEC can be carried over wholesale to the Scri-ANEC. In particular: 

• For null geodesies that come in from spatial infinity and return to spatial infinity (a > SVSM; r > 3M) the 
integrand is everywhere positive and so the ANEC holds on this entire class of geodesies. (Note that satisfaction 
of the ANEC does not arise from the trivial observation that the heat bath contributes an asymptotically 
constant and positive energy density far from the black hole — rather one has the stronger statement that the 
integrand itself is positive along the entire geodesic.) 

• The fact that the integrand is positive implies that ANEC will also be satisfied for curves with transverse 
smearing 0, or curves with arbitrary positive weighting functions |pl]] . 

• For infalling or outfalling null geodesies that reach asymptotic spatial infinity (a < SVSM; r g (2M, oo)) the 
integrand is everywhere positive and ANEC is satisfied. [Caveat: I stop the ANEC integral once it reaches the 
horizon.] 

• This may be interpreted as follows: Because Scri-NEC is satisfied, it automatically follows that Scri-ANEC is 
satisfied. 

• For the unstable circular photon orbit (a — 3^/3M; r = 3M) the integrand is everywhere positive and ANEC is 
satisfied. 

• However, for trapped null geodesies (a > 3\/3M;r £ [2M,3M]) the integrand is no longer necessarily positive, 
and this is the only case for which we will explicitly need to look at the "weighting function" ^^dA appearing in 
the ANEC integral. 



B. ANEC on trapped null geodesies? 

Observe that the "weighting function" appearing in the ANEC satisfies |6[ page 133, equations (12.60)-(12.63)] 

^^dX = goo dA = (1 - 2M/r) ( ^—r) (1 - '2M/r)dt = dt. (64) 



dXj ' V 1 - 2M/r 

Although this calculation was carried out for the Schwarzschild geometry the result remains true for an arbitrary static 
spacetime — the ANEC integral is simply the time-average along a null geodesic of the local-Lorentz NEC integrand 
where the time average is to be taken with respect to the natural static time coordinate. 

For actual calculations it is much more practical to re-express this as an integral with respect to the radial variable 
r by using 

dr 



^\ - 2M/r 

so that 



cos V' dt a/1 - 2M/r, (65) 



dr dr 

dt = = , (66) 

cosV' (l-2M/r) (1 - 2M/r) ^1 - (1 - 2Af/r)a2/r2 

Now this observation appears to weight the region near the event horizon very heavily — because of the explicit pole 
at r = 2M . However, the integrand n(r, a) has a zero at the event horizon — in the Page approximation one discovers 
an explicit factor of (1 — 2M/r) [CJ. equation |5^, while the numerical data also exhibits a first-order zero in H(r, a) 
at the horizon. Thus for trapped null geodesies one may write 

H(r,a) dr 

Ua)^2 ^ , — (67) 

' Am {l-2M/r) ^1 - (1 - 2Af/r)aV^' 

and have some confidence that the integral actually converges at the lower bound r = 2M . As a penultimate step, 
recall that calculating rniax(a) involves solving a cubic. It is more convenient to parameterize the trapped geodesic 
by calculating the impact parameter in terms of the maximum height attained by the null geodesic: 



15 



So for these trapped null geodesies one finds the ANEC integral is 

To actually evaluate this integral numerically it is useful to change variables to z = 2M/r, and zo — 2A//rmax, with 
the result that 

^.(-o) = 4M (70) 
With a little work, the square root can be seen to factorize explicitly 

/,(zo) = 4Mzo^/r^ /' ^Jf^ , ^ = ^i^, (71) 

Jz„ - 2) ^/^^ + Zl + ZZo - Z - Zq ^Jz-Zq 

This integral, though singular at the lower limit (z = zq; corresponding to r = r,„ax), is now certainly finite. While 
Mathematica has resources to deal with l/yfz singularities at the endpoints of the integration range, it reacts badly 
to such singularities when the location is chosen dynamically. That is: /^^ l/-\/z — zq is handled badly. For this 
reason the change of variables z — zq + w, while being a formal mathematical identity, leads to much better numerical 
behaviour for the integral: 

/ (zo) = 4Mzq^/1 - Zo / 7 ^^TT- — — r , —j=. 72 

Jo (zo + wy{l-w~ Zo) ^{zn + w)2 + z^ + (zq + w)zo ~ {zq + w) - zq 

This integral was determined numerically both for the Page approximation and for the Anderson-Hiscock-Samuel 
numeric data. The results are graphed in figure ^. Remember that it is meaningless to force this particular integral 
out of the range r S [2M, 3M], corresponding to zq E [2/3, 1] 




Thus we have (finally!) managed to characterize the precise class of null geodesies on which the [truncated] ANEC 
is violated. 

The critical values of zq, rmax, and a are, for the Page approximation: 

Zo = 0.90083; 
W = 2.2202 M; 

a = 7.0501 M. (73) 
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For the numerical data, one obtains: 



zo = 0.8497; 
r„,ax = 2.354 M; 

a = 6.071 M. (74) 

Note that these numbers are quahtatively reasonable and in agreement with the violations of the point-wise energy 
conditions, and the numerical investigation of the PNEC. Furthermore the ANEC violations in the numerical data 
are seen to extend out to larger distances than the ANEC violations in the Page approximation, in agreement with 
the general trend. 

Though I have quoted four-significant-digit accuracy for the numerical data one should probably not take anything 
past the second significant digit too seriously. In fact, given the vagaries of numerical analysis on singular integrals 
(and the inherent uncertainties in the Anderson-Hiscock-Samuel data) it is conceivable that the exact critical impact 
parameters might be a = 7M and a = 6M respectively. 

C. ANEC on non-geodesic curves? 

If one considers arbitrary non-geodesic curves the ANEC loses much (though not quite all) of its power: by a 
judicious choice of non-geodesic null curve one could try to remain in a region where the NEC is violated and thereby 
"trivially" violate the ANEC |l|l. 

For instance, along any circular null curve at fixed r {not a geodesic except in the case of r = 3Af) the ANEC 
integral is still proportional to p + p. Inspection of the numeric data indicates that p + p < for r < 2.298 M. (Page's 
analytic approximation giving a slightly different result, p + p < for r < 2.18994 Af.) Thus the ANEC is not satisfied 
for this particular class of non-geodesic null curves. 

XII. DISCUSSION 

Investigation of the properties of the averaged null energy condition is of considerable interest to diverse applications 
in semiclassical quantum gravity. It is now abundantly clear that, in the test-field limit, semiclassical quantum fields do 
noi generally satisfy the ANEC: Indeed ANEC violations are related to the existence of a non-zero scale anomaly []6|,p0|. 
Even if the scale anomaly vanishes, this does not necessarily imply that the ANEC is satisfied: one has to do a case 
by case analysis. As an example, this paper investigates the situation in Schwarzschild spacetime. 

The analysis presented herein is somewhat of an attempt to crack a walnut with a sledgehammer in the sense that 
it is a collection of rather general techniques applied to a rather particular problem — it is certainly true that this 
type of analysis can now be carried forward to other geometries, other quantum fields, and other vacuum states by 
by straightforward but tedious computation. 

For the Schwarzschild geometry (with a conformally coupled scalar field in the Hartle-Hawking vacuum) the results 
may be expressed thusly: 

• Point-wise energy conditions — 

— Inside the event horizon, with suitable caveats regarding the applicability of the Page approximation, 
almost any energy condition you can think of will be violated. 

— Between the event horizon (r = 2M) and the unstable photon orbit (r = 3A/) many of the energy conditions 
are violated, in a series of onion-like layers. 

— Outside the unstable photon orbit (r — 3M) all energy conditions are satisfied. 

— All null geodesies that reach asymptotic infinity are well-behaved: If you look along the null geodesic you 
never see NEC violations. 

• Averaged null energy condition — 

— If your null geodesic ever reaches infinity, the [truncated] ANEC is definitely satisfied, and is satisfied for 
a non-trivial reason: the integrand is strictly positive all the way from the event horizon to null infinity. 

— Between the event horizon (r — 2M) and the unstable photon orbit [r = 3M) some of the trapped null 
geodesies violate the [truncated] ANEC. 
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— If you are willing to look at non-geodesic null curves even more violations of the ANEC can be found. 
It should be possible to generalize the observations of this paper. For instance: 

1. It would be very nice to have an analytic understanding of the precise role played by the unstable photon 
orbit — the numerical evidence is suggestive but not definitive. 

2. For that matter it would be nice to know if the relevance of the unstable photon orbit generalizes to other 
geometries. 

3. It would be nice to go beyond the numerics; to develop some exact analytic arguments that go beyond the Page 
approximation. 

4. Generalizations to the Boulware vacuum will be presented in a companion paper. 

5. Generalizations to the Unruh vacuum, other quantum fields, non-conformal couplings, particle masses, and the 
Reissner-Nordstrom geometry will be straightforward if tedious. 

6. The new energy conditions I introduce, Scri-NEC and Scri-ANEC, are interesting in that they focus attention 
on null infinity. And null infinity is where all the interesting details arise in the Friedman-Schleich-Witt |2|] 
topological censorship theorem, and the Penrose-Sorkin-Woolgar version of the positive mass theorem |2^] . I 
suspect that with a little more work suitable generalizations of these theorems can be constructed in terms of 
the Scri-ANEC. 

7. Finally I should point out that even though the various energy conditions are violated in many regions, this 
does not give one a completely free hand to design spacetime geometries to taste: It seems quite likely that 
the "quantum inequality" approach of Ford and Roman p9| , ^ will allow us to place constraints on spacetime 
geometries even if all the usual types of energy condition fail. 
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